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ABSTRACT 

In this paper we develop an analytical theory for the interaction between waves and currents 
induced by breaking waves on time-scales longer than the individual waves. We employed the wave- 
averaging procedure that is commonly used in the literature. The near-shore zone is often characterized 
by the presence of breaking waves. Hence we develop equations to be used outside the surf zone, based 
on small-amplitude wave theory, and another set of equations to be used inside the surf zone, based on 
an empirical representation of breaking waves. Suitable matching conditions are applied at the 
boundary between the offshore shoaling zone and the near-shore surf zone. Essentially we derive a 
model for the interaction between waves and currents. Both sets of equation are obtained by averaging 
the basic equations over the wave phase. Thus the analytical solution constructed is a free vortex 
defined in both shoaling and surf zones. The surf zone solution is perturbed by a longshore component 
of the current. Thus the presence of the rip current cell combined with the longshore modulation in the 
wave forcing can drive longshore currents along the beach. The outcome, for our set of typical beach 
profile, is a description of rip currents. 

KEYWORDS'. Wave-current interactions, surf zone, shoaling zone, matching conditions, wave- 
averaging, rip currents, radiation stress. 



I. INTRODUCTION 

The action of shoaling waves, and wave breaking in the surf zone, in generating a wave-generated 
mean sea-level is well-known and has been extensively studied, see for instance the monographs of Mei (1983) 
and Svendsen (2006). The simplest model is obtained by averaging the oscillatory wave field over the wave 
phase to obtain a set of equations describing the evolution of the mean fields in the shoaling zone based on 
small-amplitude wave theory and then combining these with averaged mass and momentum equations in the 
surf zone, where empirical formulae are used for the breaking waves. These lead to a prediction of steady set- 
down in the shoaling zone, and a set-up in the surf zone. This agrees quite well with experiments and 
observations, see Bowen et al (1968) for instance. However, these models assume that the sea bottom is rigid, 
and ignore the possible effects of sand transport by the wave currents, and the wave -generated mean currents. 
Hydrodynamic flow regimes where the mean currents essentially form one or more circulation cells are known 
as rip currents. These form due to forcing by longshore variability in the incident wave field, or the effect of 
longshore variability in the bottom topography (Kennedy 2003, 2005 ,Yu & Slinn 2003, Yu 2006 and others). 
They are often associated with significant bottom sediment transport, and are dangerous features on many surf 
beaches (Lascody 1998 & Kennedy 2005). 

There is a vast literature on rip current due to wave-current interactions, see the recent works by 
(Horikawa 1978 , Damgaard et al. 2002, Ozkan-Haller & Kirby 2003 ,Yu & Slinn 2003 , Yu 2006, Falques, 
Calvete & Monototo 1998a and Falques et al 1999b, Zhang et al 2004 and others) and the references therein. 
Our purpose in this paper is to exploit the fully developed but under-utilised wave-current interaction theory in 
the nearshore.In section 2 we record the usual wave-averaged mean field equations that are commonly used in 
the literature. In section 3 we introduce a description of the rip current formation and examine the consequences 
for both shoaling and surf zones. Then in section 4 we employ section 3 to a choice of linear depth profile. We 
conclude with a discussion in section 5. 

II. FORMULATION 

2.1 Wave field 

In this section we recall the wave-averaged mean flow and wave action equations that are commonly 
used to describe the near-shore circulation (see Mei 1983 or Svendsen 2006 for instance). We suppose that the 
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depth and the mean flow are slowly-varying compared to the waves. Then we define a wave-phase averaging 
operator < / >= / , so that we can express all quantities as a mean field and a wave perturbation, denoted by a 
"tilde" overbar. For instance, 

where C, is the free surface elevation above the undisturbed depth h = h(x) . Then outside the surf zone, the 
representation for slowly- varying, small-amplitude waves is, in standard notation, 

C,(x,i) : a cos 6. (2) 
Here a = a(x.t) is the wave amplitude and 6 = 0(x, t) is the phase, such that the wavenumber k , frequency 
Q are given by 

k = (k,l) = V0 Q = -6>. (3) 

The local dispersion relation is 

Q. = co + k.U, a> 2 = grc tanhxf/ (4) 
where k 2 = k 2 +l 2 . 

Here U (x,t) is the slowly- varying depth-averaged mean current (see below), and H(x,t) = h(x) + £(x,t) . 
To leading order, the horizontal and vertical components of the wave velocity field are respectively 
~ . k coshfc(z + h) smh/c(z + h) . 

u .—ail cost/, w .ail sine'. (5) 

k sinh kH sinh kH 

Importantly, note that we have ignored here any reflected wave field, which is assumed to be very weak when 

the bottom topography is slowly varying. 

The basic equations governing the wave field is then the kinematic equation for conservation of waves 

k t +Vco = 0, (6) 

which is obtained from (3) by cross-differentiation, the local dispersion relation (4), and the wave action 
equation for the wave amplitude 

A+V-(c g A) = 0. (7) 

Here A = ElcO, where E = ga 2 /2 is the wave energy per unit mass, and 
C g =V k CO = U + C g k/K,(c g = dcoldK) is the group velocity. Using the dispersion relation (4) in (6) we 
get 

k t +c g -Vk = -V e P, (8) 
Q,+c e -VQ = Q t , r (9) 
Here the subscript "e" denotes the explicit derivative of Q(k,X, t) with respect to either X or t , when the 
wavenumber k is held fixed. In this water case these explicit derivatives arise through the dependence of Q on 
the mean height H and the mean current U . 

2.2 Mean fields 

The equations governing the mean fields are obtained by averaging the depth -integrated Euler 
equations over the wave phase. Thus the averaged equation for conservation of mass is 

H t +V(HU) = 0. (10) 

Here H = h + £ where h = h(x) is the time-independent undisturbed depth. For the velocity field we 
proceed in a slightly different way, that is we define 

U = U+U, (11) 
where we define U so that the mean momentum density is given by 

M = HU =< udz >, (12) 

J-h 

But now we need to note that u does not necessarily have zero mean, and that U and U are not necessarily 
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the same. Indeed, from (11) and (12) we get that 

rC ■ 

u = U+ <u >, and < \ u dz >= 0. 

J-h 

But U =U + 0(a 2 ) , so that < U > is 0(a 2 ) and it follows that, correct to second order in wave amplitude, 

E k 

M=Hu + M w , where M w = -H < u >=<£u(x,0,t) >= . (13) 

C K 

The term M in (13) is called the wave momentum. 

Next, averaging the depth-integrated horizontal momentum equation yields (Mei 1983) 

(HU ). + V.(HUU )=- V. < f uu + pldz > + <p(z = -h) > Vh. 

J-h 

Next an estimate of the bottom pressure term is made by averaging the vertical momentum equation to get 

- rC K 

< P(z = -h)>-g(£ +/i)=V. < wudz> + <\ wdz> r 

J-h J-h 

For slowly-varying small-amplitude waves, the integral terms on the right-hand side may be neglected, and so 
< p(z = —h) >~g(£ + h) . Using this in the averaged horizontal momentum equation, and replacing the 
pressure p with the dynamic pressure q= p + (z~^) yields 

HU t +V.(HUU)=-V.S -gHV^ (14) 

K g ~? 

where S=<\ [uu + qI]dz>-<—£ > I. (15) 

J-h 2 

Here S is the radiation stress tensor. In the absence of any background current, so that U is 0(a 2 ) , we may 
use the linearized expressions (2, 5) to find that 

S*c g k- + E[^-^I. (16) 
8 co c 2 

where the phase speed C = CoIk , correct to second order in the wave amplitude. 
2.3 Shoaling zone 

These equations hold in the shoaling zone outside the surf zone (defined below). In summary, the 
equations to be solved are that for the conservation of waves (6) combined with the dispersion relation (4), the 
wave action equation (7), the averaged equation for conservation of mass (10) and the averaged equation for 
conservation of horizontal momentum (14), where the radiation stress tensor is given by (16). In this shoaling 
zone, we assume that wave amplitudes are small, and that there is no background current. Then all mean 

quantities are 0(a 2 ), and in particular we can systematically replace H with h throughout these 
equations. Next we shall suppose that h = h(x) depends only on the offshore co-ordinate X > , where the 
undisturbed shoreline is at X = defined by h(0) = . Further, in the near-shore region, including all of the 
surf zone, we shall assume that h > . Then we seek steady solutions of the equation set in which all variables 
are independent of the time t , and are also independent of the transverse variable y . It then follows from the 
mean mass equation (10) that HU is constant, and since H = at the shoreline, it follows that we can set 
U = everywhere. Then in the dispersion relation (4) Q = CD . From the equation for conservation of waves 
(6) we see that the frequency CO and the transverse wavenumber / are constants, and the the offshore 
wavenumber k is then determined from the dispersion relation (4). As is well-known, it then follows that as 
H — > , I k I— > go , that is the waves refract towards the onshore direction, where we assume that the waves 

are propagating towards the shoreline so that k < 0. The wave action equation (7) reduces to Ec is constant. 
Near the shore, we can assume that the shallow water approximation holds and then C g ~ —(gh) l/2 , so that 

an * a /J , (17) 
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where a is the wave amplitude at a location offshore where h = h Q . The surf zone X < X b , h < h h = h(x b ) 
can then be defined by the criterion that h b is that depth where alh = A cr , defining an empirical breaking 
condition. A suitable value is A cr = 0.44 , see Mei (1983) or Svendsen (2006). 

The last step is to find the wave set-up C, from the mean momentum equation (14), which here 

becomes 

S x + gHC x =0„ (18) 

where S = ^E CO s 2 <p + (—--)E. 

c c 2 

Here (/) is the angle between the wave direction and the onshore direction, and S is the " XX " component of 
the tensor S . As h —> , <f) ~ 0, C g ~ C, S ~ 3E/2 , and we recover the well-known result of a wave set- 
down in the shoaling zone 

2 2 i 1/2 

Z=- a ^ = -^- (19) 

Ah Ah 

Here we have assumed that C, is zero far offshore. Note that the first expression for Q does not need the use 
of the shallow water approximation, as shown by Longuet-Higgins and Stewart (1962). 

2.4 Surf zone 

In the surf zone < X < X b ,0 < h < h b , we make the usual assumption (see Mei (1983) for instance) 
that the breaking wave height 2a is proportional to the total depth H , so that 

2a = yH, or E = ^- , (20) 

8 

Here the constant y is determined empirically, and a typical value is y = 0.88.. To determine the mean 
height H = h + £, we again use the mean momentum equation (14), but now assume that 
S = 3E/2 = TgH 2 l2 where T = 3y 2 /S , so that 

THH+H(H-h)=0, so that H = H h + ~\ (21) 

(l + D 

where the constant H b = h b + Q h is determined by requiring continuity of the total mean height at X = X b . 
Note that using (19) 

H b = h b -a 2 X 2 IMl 12 , 

and since H b must be positive, there is a restriction on either the deep-water wave amplitude a or on the 
breaker depth h b , 

allAhl<hflhf. (22) 
Note that the expression (21) is valid for any depth h{x) , although in the literature it is often derived only for a 
linear depth profile h = OCX . 

We are now in a position to determine the displaced shoreline X = X s , defined by the condition that 

H = . That is, if h s = h(x s ) then H = (h - h 5 )/l + T , where 

h s =-Th b -(l + T)£ b , (23) 
Note that to use the expression (23) it may be necessary to extend the definition h(x) into X < . For instance 
for a linear beach, h = OCX , this is straightforward, but for a quadratic beach profile, h = /3x 2 , the extension 
for negative X should be h = —fix say. Note that from (19), 
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a 2 h m 



Ah b 



3/2 



and, on combining this with the condition (22), it follows that the shoreline recedes (advances), that is 
h, < 0(> 0) when 



2 -p 7,5/2 

V < . F k b 



Ah 2 1 + T h*' 2 ' 

, 5/2 2 

r h, a A 



or —^ < ^ <h 5 b p -/h*> 2 . (24) 

l+r/jf 4h 2 b 

Curiously, this anomalous result does not seem to have been noticed previously [5] and [5]. Since there is an 
expectation that the shoreline should advance (see Dean and Dalrymple (2002) for instance), essentially it states 
that the present model is only valid for sufficiently small waves far offshore, defined by the first inequality in 
(24), which slightly refines the constraint (22). 

III. A GENERAL DESCRIPTION OF THE RIP CURRENT FORMATION 

Here we consider a steady-state model driven by an incident wave field which has an imposed 
longshore variability. The wave field satisfies equation (7) which in the present steady-state case reduces to 

ddx(Ec g cos 9) + ddy(Ec g sin 9) = 0. (25) 
Here we again assume that h = h(x) and that consequently the frequency CO and the longshore wavenumber 
/ are constants, while the onshore wavenumber K is then determined from equation (4). We have the wave 
energy E of the form 

E = E Q (x) cos Ky + F Q (x) sin Ky + G (x), (26) 
where the longshore period ItiIK is imposed. These equations in the shoaling zone yields 

(E c g cos 0) x + KF Q c g sin 9 = (27) 

(F c g cos 9) x - KE c g s in 9 = (28) 

(G c g cos 9) x =0. (29) 

on collecting terms in COsAj" , sinA/y and the constant term, which form three equations for E , F Q and 

G . Equation (29) easily yields that G C COS 9 = constant. In shallow water, we may approximate by 

putting C g « gh m and COSr? ~ 1 , so that then G ~constant/h U2 . For the remaining equations we can use 

Snell's law, sin 91c = sin 9 b lc b = CX b (the constant value, here evaluated at the breaker line), and the shallow- 
water approximation to get that 

{(E Q c) x )/c} x + K 2 a 2 b E Q c 2 = 0, (30) 

while although F satisfies the same equation, once E has been found, then F is given by either (27) or 

(28). In practice, Kc « 1 and so approximately we can assume that (E Q , F Q )c ~COnstant , the usual 

shallow- water expressions. Note that here C » -\] gh . In the surf zone, the expressions E Q (x), F Q (x), G Q (x) is 
determined empirically. 

Once the expression (26) has been determined, we may then substitute into the expressions (27,28 & 
29) to obtain the radiation stress fields. Our aim here then is to describe how steady-state rip currents are forced 
by this longshore modulation of the incident wave field, especially in the surf zone. 

The forced two-dimensional shallow water equations that we use here are characteristic of many 
nearshore studies (Horikawa 1978 , Damgaard et al. 2002, Ozkan-Haller & Kirby 2003 ,Yu & Slinn 2003 , Yu 
2006, Falques, Calvete & Monototo 1998a and Falques et al 1999b, Zhang et al 2004 and others). Then, 
omitting the overbars as before, then equations (??) in the present steady-state case reduce to 

H[UdUdx + VdUdy] = -gHd<@x-[T x ], (31) 
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H[UdVdx + VdVdy] = -gHd^dy - [r y ], (32) 
where the stress terms are defined; 

t x = dS n 8x + 8S 12 dy andr y = 8S 21 dx + 8S 22 8y. (33) 

Next we observe that equation (10) can be solved using a transport stream function y/(x, y) , that is 

U = — —dy/dy and = — dy/dx, (34) 
H H 

Next, eliminating the pressure, we get the mean vorticity equation 

V'^y "^' ( H )r = [ 7/ ]y ~ [ il ]x (35) 

where Q is define as 

n = V x -U y =(^) x + (^) y . (36) 

We shall solve this equation (35) in the shoaling zone X > X b and in the surf zone X < X b , where as before 

X = x b is the fixed breaker line. It will turn out that the wave forcing occurs only in the surf zone, but 
continuity implies that the currents generated in the surf zone must be continued into the shoaling zone. 

3.1 Shoaling zone 

In X> X b we shall assume that H ~ h as £ is 0(a 2 ) . Then we shall use the expressions [27 ,28] 
to evaluate the radiation stress tensor. For simplicity, we shall also use the shallow-water approximation that 
C g ~C « ^jgh , and so we get that 

S u =E( co i9 + ^),S n = S 2l =£sin6>cos6>,S 22 =£(sin 2 6> + ^) (37) 

These expressions are in principal known at this stage, and so we can proceed to evaluate the forcing term on 
the right-hand side of (35). To assist with this we recall Snell's law 

sin # = yfhJh. sin h 



where h b and 6 b are the water depth and incidence angle at the breaker-line. Now the energy equation (25) has 
the approximate form 

(£ccos G) x + (£csin 0) y = 0, 
and using Snell's law, this can be written as 

(£cos 2 6>) v +(£sin#cos6>) v + £— = 0, 

c 

and so t=—E-E—. 

2 c 

We can also deduce from (25) that 

(Esm0cos0) x +(Esm0) y =O, 

and so r„ = —E v . 

- 2 3 

We can now evaluate the right-hand side of (35), and find that its identically zero, 

[— ] -[— ] =o. 

l h y h h 

Thus in the shoaling zone there is no wave forcing in the mean vorticity equation, although of course there will 
be a mean pressure gradient. However, this does not concern us since here our aim is to find only the flow field. 
Note that the result that there is no wave forcing in the vorticity equation does not need the specific form (26), 
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and is based solely on the steady-state wave energy equation (25). The specific form (26) is only used in the surf 
zone. 

With no forcing term, the vorticity equation (35) can be solved in the compact form, noting that we again 
approximate H with h , 

^ = F(y/). (38) 
h 

But here F(y/) = from the boundary conditions in the deep water as X — > oo , where the flow field is zero. 
Thus our rip current model has zero vorticity in the shoaling zone. It follows that we must solve the equation 

n = (-y/ x ) x HT¥y) y =Q, 09) 

h h 

in X > X b . Since h = hix) we can seek solutions in the separated form 

y/ = X{x)Y{y) (40) 
with the outcome that 

(4^-^ = 0, Y yy + K 2 Y = 0. (41) 
h n 

We note the separation constant K 2 = 2tz/L must not be zero, and is in fact chosen to be consistent with the 
modulation wavenumber of the wave forcing. Without loss of generality, we can choose 

Y = sinKy. (42) 

For each specific choice of h{x) we must then solve for X(x) in X > X b , with the boundary condition that 
X — > as X — > co . We shall give details in the following subsections. Otherwise we complete the solution by 
solving the system (35) in the surf zone, and matching the solutions at the breakerline, X = X b where the 
streamfunction y/ must be continuous, and in order to have a continuous velocity field we must also have that 
l// x is continuous. 
3.2 Surf zone 

To make sense of wave forcing, we assume that the expression (26) holds in this region. The functions 
E (x), F Q (x), G (x) are then determined empirically. To determine the wave forcing term in the mean 

vorticity equation (35) we shall assume that 9 = 9 h « 1 so that, on using (??) and (37) we get that 

T = —E , T = — E . 

x 2 «• v 2 V 

Then (35) now becomes, where we again approximate H with h(x) , 

E xy , E y K = (h m E v \ 
h 2h 2 h 3 

where here Q = QJh is the potential vorticity. Since the wave forcing is given by (26), that is 

E = E (x)cosKy + F Q (x) sinKy + G (x), (44) 

we observe that the unmodulated term G () (x) plays no role here at all, although of course it will contribute to 

the wave setup. In order to match at X = X b with the expression (42) for the streamfunction in the shoaling 
zone, we should try for a solution of (43) of the form 

y/ = F(x)smKy + G(x), in x<x b . (45) 

The matching conditions for the streamfunction and velocity field at the breakerline X = X b require that 

F(x b ) = X(x b ),F x (x b ) = X x (x b ),G(x = x h ) = 0, G x (x„ ) = 0. 
The expression (45) yields 

Q = F sinKy + G (46) 

where F and G are differential operators where they are defined as; 



8 y/dxdCldy - d y/dydQdx = ^ + = (43) 
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F K^F 

F = (~f) x T- (47) 

h h 

G = Z X ,Z = ^ (48) 
n 

From equation (43) we get a set of three equations that are used to determine the rip-current flow field in the 
surf zone. These are namely; 

^T- F( T^ =0 ' (49) 
h h 

1/2 



(G--F(-) x )= K iy » J -\ (50) 
(^^ = 0. (51) 



F(CG-(-) x )= K (53) 



h h h 
h 

Equation (51) gives E '. 1/h 112 , which is an unacceptable singularity as h — >0 . Hence we must infer that in 
the surf zone at least, E = . The first of the three equations, that is (3.2a) suggests that 
F 

— = CF whereCisaconstant, (52) 
h 

and the second (3.2b) yields that 

r G_ )= (^F ) X 
h >x> h 

The boundary conditions at X = where h = are that both mass transport fields U,V should vanish, that 
is from (34) Xj/ = constant and y/ Jh = , which implies that 

F = F K =0,G = constant, — - = 0, at x = 0. (54) 

h 

As above there are also the matching conditions for both F and G separately at the breakerline, that is for F 
we have that 

F(x„) X(x t ) 

where we note that here the right-hand side is a known quantity, depending only on K and X b . Next we see 
that equation (52) reduces to 

F K 2 F 

{-f) x — = ChF (55) 

h n 

Together with the boundary conditions at X = 0, X = X h this is essentially an eigenvalue problem for F(x) 

with eigenvalue C . In general it is solved approximately since we shall assume that Kx b « 1 . Once F(x) is 

known we can solve (53), together with the appropriate boundary conditions to get G(x) to complete the 
solution. 

Note that the amplitude of F(x) is an arbitrary constant in this solution, and so we can fix it by 

specifying its value at X = X b say. Indeed the solution we have constructed is essentially a free vortex defined 

by X(x)sin^y in the shoaling zone X> X b , and F(x)sklKy in the surf zone X< X b , perturbed by a 

longshore component G(x) in the surf zone. Note that in the presence of the wave forcing, both F, G are 
non-zero, see (53). It is significant that unlike the longshore currents considered in the basic state which depend 
on an ad hoc fractional parametrization, the presence of the rip current cell combined with the longshore 
modulation in the wave forcing can drive a longshore current. 
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IV. AN APPLICATION TO LINEAR DEPTH PROFILE 

To make sense of the rigorous mathematical derivations and formulations we then suppose that 
h = OCX . Then immediately the solutions of equation (41) in the shoaling zone is 

X(x) = C x xK, (Kx) + C xl, (Kx). 

Since /j becomes unbounded at X = oo , the boundary condition X — > as X — > oo gives C = , and so 
we get that 

X(x) = C 1 xK l (Kx),x> x b . (56) 
The behavior of the Bessel functions depict the rip-current character which decays exponentially away 
from the breakeiiine, on the scale K 1 . It is interesting to note here that the bottom slope CC does not appear in 
this solution at all. In the surf zone, the boundary condition (54) implies that as X — > , F(x) ~ AqX 2 . We 
shall make the approximation that the shoaling zone has a small cross-shore width, in which case it is easily 
seen that the right-hand side of (55) is much smaller that the left-hand side, and is 0(x 3 ) . Approximating the 
right-hand side accordingly we get that 

F xx --F x -K 2 F = \Ax\ (57) 
x 

where A = Ccc 2 . This is useful, as the solutions of the homogeneous equation on the left-hand side are known, 
namely j = xI { (Kx) and j b = xK x (Kx) . By the variation of parameters one gets 

F(x) = A(x) j a + B(x) j b + C, j a + C 2 j b , (58) 

noting that W(j a , j b ) = —X , where 

A(x)= [ j b A Ax J 'dx, 5(x) = -[ j^Ax^dx. (59) 

JO JO 

The boundary condition at x = shows that C 2 = 0, and the normalization of F as X — > implies that 
kC { = 2A Q . It remains to apply the boundary condition at X = X b which then yields A (that is C). From 
equations (58), (59) one gets 

F«V 2 (1 + ^)- 

Note that to this order, F is independent of K where A scales as x b so this small X b approximation also 
requires that the constant A also be very small. Similarly the right-hand side may be approximated by 

K 2 x;(r +MKx h /2)) 
where y Q is Euler's constant. This implies that to leading order 

Ax b a -8. (60) 
But this is a bit too simple as can be seen from the consequent expression for 



FvAvxHl-^—), (61) 



x 4 



3x b 

so that actually F « at X = X b , which is too simple. The higher-order terms in A and B can be found, 
and then we get that 

F = \{x 2 --x 4 + — x 6 + A[— + — - — K 2 x s ]). 
* 8 192 24 128 72 

We define — Axt + 1 = £ « 1 and keeping the required next order terms we get that 
8 

e = ^K 2 x 2 b [(r +MKx b )) + ^l 
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Figure 1: Plot of i^xY^, against xlx b , where A is arbitrary as given by equation (61). The figure shows 
that there is only curve at leading order. We note that there is no dependence on the slope CC and only a weak 
dependence on K .We observe that F(x) reaches a maximum value of F(x)/A Q ~ 0.7. Beyond this point 

F(x) may probably decrease monotonically. 
Next for the y -independent component, we need to solve for G(x) from(53). As above, we 
approximate F = A^x 2 , and we use the empirical expression F = y 2 h 2 /8 . This yields 

Z xx --Z x -a 2 x 2 CG = (62) 
x \6AqX 

There is a singularity at x = , which can be analyzed by setting U = X 2 so that (62) becomes 

4Z„„-^Z = -^4 (63) 

where we have put X = —A, > and the full solution is 

„ 5 ga 2 y 2 x h , x 1 . r V2x 2 nN 

Z= ' ( r^sm[^-]) (64) 

12^ x 6 sinV2 x fc 

Finally the complete the stream function G(x) is found by integrating G v = hZ subject to the boundary 

condition that G = at X = X b . Thus we get that, for < X < X b , 

G(,x) = — — + ^ ^cos[^^ r (65) 

124, 3x fc 3 2V2sinV2 x 2 3 2V2tanV2 

Note in particular that G(0) # and is the net mean longshore mass transport in the rip current system. 



Figure 2: Depicts the plots of normalized Z(x) and G(x) given by equations 64 and 65 where each is 

normalized by 5gCC 2 y 2 X b l\1A Q and 5 g CC 2, y 2 1 \1 respectively with A Q arbitrary. Observe that in the 
figure [2] there is a small region of reversed flow near the breaker line. 

The combined expressions (56, 61, 65) complete the solution, where we recall that the constant C is 
given by (60) (since X = Coc 2 '), or their respective higher-order corrections. Note that the amplitude of F(x) 
at X = x b is given by 

F(x b ) = C l x b K l (Kx b ). 
On using the approximation Kx b « 1 , and the approximate expression (61), this reduces to 



F(x h ) ■■ 



2A)X b _ C[ 



3 K 

The rip-current system contains a free parameter A Q or its equivalent. We choose to define this free parameter 
to be the value of F{x b ) and normalize the full solution by this value. Thus we get from (40, 42) in X > X b , 
and (45, 61, 65) in X < X h that the normalized streamfunction XjJ n is given by 
X(x) 



Wn 



X(x h ) 



sin(Ky), for x>x b , 



(66) 
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Fix) G(x) 

^«=^H sini ^ + 7? 7^' f° r < x < x b- (67) 
F(x b ) G(0) 

Here R = G(0)/ F(x b ) is a free parameter. Thus a larger (smaller) R decreases the circulation of the rip- 
current system vis-a-vis that of the longshore current component. From (61, 65) we find that here 

R = 5gaYx 1 1 

8^ 3 2V2tanV2 ' 
Note that with all other parameters fixed, a larger (smaller) slope Ct increases (decreases) R . In order to 
estimate typical values for R we note that from (61) the longshore velocity field in the " SVCl(Ky) "-component 

scales as V c = \J Ct , while the longshore component then scales with RV C . Taking account of the actual 

numerical values in the expressions given above, we find that a suitable values are R » —0.1. Plots of y/ n are 

shown in figures 3, ??, 4 for R = —0.02,-0. 1,-0.2,— 0.5 —2 respectively, with Kx b = 0.2 (noting that our 

present theory requires that Kx b is small). 

Figure 3: Plot of the rip current streamlines for a linear depth profile, given by equation (67) where F(x) and 
G(x) are equations (61) and (65) respectively for R = —0.02 in the left panel and R = —0. 1 in the right 

panel 

Figure 4: As for figure 3 but R = —0.5 in the left panel and R = —0.2 in the right panel 

From the plots in figures 3, ?? we see that as I R I increases, the core of the rip current circulation 
moves from the shoaling zone towards the surf zone. The reason for this is that the solution we have constructed 
is essentially a free vortex defined by X(x)sin Ky in the shoaling zone X > X b , and F(x)sh\ Ky in the surf 

zone X < X b , perturbed by a longshore component G(x) in the surf zone. Since R < , this longshore 
component opposes the vortex flow in the cell < Ky < n but is in sympathy for n < Ky < In . This has 
the effect of moving the vortex cell further offshore in the sector n < Ky < In relative to the sector 
< Ky < n . Note that I R I increases as the wave forcing increases, or as the slope Ct increases, or as the 
depth CCX b at the breaker line increases. 

V. CONCLUSION 

We described qualitative solutions for rip currents which are essentially free vortices in both zones. 
The free vortex in the surf zone is perturbed by a longshore modulation in the wave forcing. Rip current cell 
combining with the longshore modulation in the wave forcing can drive longshore currents along the beach. 
Thus the dynamics of the shoaling zone is only dependent on the state-state wave energy equation.The wave 
forcing in the surf zone sets the wave activities different from those of the shoaling zone. To determine wave 
forcing in the mean vorticity equation we assume that the wave angle becomes smaller. We also note here that 
the component of the radiation stress in the y momentum remains unchanged across the entire flow domain. 

This shows that it is only the X component of the radiation stress that play a leading role in the wave forcing. 
However, wave forcing encountered in the surf zone has an unmodulated term that does not play a role in the 
vorticity equation but only contribute to wave setup. To ensure continuity of the streamfunctions in the shoaling 
zone we match the solution at the breakerline by a matching condition with appropriate boundary conditions. 
Thus the rip currents solution in the surf zone is provided by the terms in the matching condition. The terms in 
the matching condition has a cross-shore width and a modulated longshore component. The cross-shore width 
was determined by the application of perturbation method and variation of parameter. It would to interesting to 
examine the effect of friction on the rip currents. 
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